Universal Neel Temperature in Three-Dimensional Quantum Antiferromagnets 



(N 

O 

(N 



u 

I 

a 
o 

in 



Songbo Jin and Anders W. Sandvik 
Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA 

(Dated: February 16, 2012) 

We study three-dimensional dimerized 5 = 1/2 Heisenberg antiferromagnets, using quantum 
Monte Carlo simulations of systems with three different dimerization patterns. We propose a way 
to relate the Neel temperature Tjv to the staggered moment ms of the ground state. Mean-field 
arguments suggest Tjv oc rus close to a quantum-critical point. We find an almost perfect universality 
(including the prefactor) if Tat is normalized by a proper lattice-scale energy. We show that the 
temperature T* at which the magnetic susceptibility has a maximum is a good choise, i.e., Tn /T* 
versus rris is a universal function (also beyond the linear regime). These results are useful for 
analyzing experiments on systems where the spin couplings are not known precisely, e.g., TlCuCls. 

PACS numbers: 75.10.-b, 75.40.Cx, 75.10.Jm, 75.40.Mg 



Quantum fluctuation can drive continuous phase tran- 
sitions between different kinds of ground states of many- 
body systems. While transitions taking place at tem- 
perature T > are controlled by thermal fluctuations, 
quantum fluctuations also play a role here. Quantum- 
critical scaling can often be observed throughout a wide 
region (the quantum-critical "fan") extending out from 
the quantum-critical point (pc, T ~ Q) into the plane 
{g,T > 0) [ll-lll, where g is the parameter tuning the 
strength of the quantum fluctuations. In addition, the 
quantum fluctuations of course also strongly affect the 
critical temperature Tc, because Tc —> as g — >■ ^c- One 
can regard the quantum fluctuations as reducing the or- 
der at low temperature (T <C Tc), with the thermal fluc- 
tuations eventually destroying it as T — > Tc, but precisely 
how the two kinds of fluctuations act in conjunction with 
each other to govern Tc is not known in general. 

We will here discuss manifestations of the interplay of 
quantum and thermal fluctuations for < T < Tn in 
three-dimensional (3D) S = 1/2 quantum antiferromag- 
nets with Heisenberg interactions. In these systems one 
can vary the critical Neel-ordering temperature, T/v, and 
ultimately achieve a quantum phase transition (Tn — s- 0), 
by considering dimerized couplings, such that each spin 
belongs exactly to one dimer and the intra- and inter- 
dimer couplings are different. The Hamiltonian for such 
models can be generically written as 



J2 E 



(1) 



where {i,j)a denotes a pair of spins coupled at strength 
Ja, with a = 1 and a = 2 corresponding to inter- and 
intra-dimer bonds, respectively. Three examples of such 
dimerized 3D lattices are shown in Fig. [TJ In (a) and (b) 
the spins form simple cubic lattices, and each nearest- 
neighbor site pair is coupled either by Ji or J2. In (c) 
two different cubes each have all Ji couplings, and pairs 
of spins in different cubes form the J2-coupled dimers. 
We will use the ratio g = J2/ Ji as the tuning parame- 
ter. When g « 1 the system is Neel-ordered at T = 
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FIG. 1: (Color online) Dimerized 3D lattices; (a) columnar 
dimers, (b) staggered dimers, and (c) double cube. For a 
system of length L, the number of spins is A'" = L'^ in (a) and 
(b), and A*' = 21/^ in (c). The two different coupling strengths 
Ji and J2 are indicated by thin (black dashed and solid) and 
thick (red) lines, respectively. 

and when g ^ cx) it decouples to form a set of inde- 
pendent dimers, with the ground state becoming a triv- 
ial quantum paramagnet with a singlet-product ground 
state. The system for any Ji > and J2 > is accessi- 
ble to unbiased numerical studies with efficient quantum 
Monte Carlo (QMC) methods with loop updates [l-Q- 

Analogous dimerized Heisenberg models have been 
studied extensively with QMC in two dimensions, where 
there is order only at T = (for g < gc) and the nature of 
the quantum-critical point and its associated scaling fan 
has been the main focus of interest 0, § . Some simula- 
tions have also been previously carried out for 3D dimer- 
ized models . Here we report calculations uncover- 
ing universal aspects of the ordering temperature, from 
systems close to the quantum-critical point to deep in- 
side the Neel phase. We develop a scaling procedure of 
direct relevance to experiments. Our results also provide 
new insights into the relevant energy scales present in the 
3D Neel state and demonstrate an effective decoupling of 
thermal and quantum fluctuations. 

Experimental issues. — The best experimental realiza- 
tion so far of a dimerized system with a quan tum phase 
transition is TlCuCla under pressure 12h14|. The spin 
dimers here form on pairs of Cu atoms that can clearly 
be identified as the most strongly coupled neighbors. 
The inter-dimer couplings are, however, more compli- 
cated than in the simple nearest-neighbor Hamiltonian 
P]) . There are several significant exchange constants but 
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their exact values are not known (although they have 
been estimated based on approximate calculations of the 
magnon dispersion, which can be compared with experi- 



ments [12|, |15[). The dimers nevertheless form a 3D net- 



work, and one can expect the same ground state phases 
and phase transitions as with the simplified Hamilto- 
nian ([1]). Under ambient pressure, TlCuCls exhibits no 
magnetic order, but beyond a critical pressure antifer- 
romagnetic order emerges continuously. The interpre- 
tation of this is that one or several of the inter-dimer 
couplings become strong enough for Neel order to form. 
The observed longitudinal and transversal excitation en- 
ergies agree well with predictions based on 0(3) symme- 
try breaking and Goldstone modes 15, lif . 

The fact that the microscopic spin-spin couplings in 
TlCuCla, and how they depend on pressure, are not 
known accurately is a complication when comparing ex- 
perimental results with calculations for a specific model 
Hamiltonian. In this situation it is useful to make com- 
parisons that do not require any explicit knowledge of the 
couplings. Here we will investigate how the Neel temper- 
ature is related to the staggered magnetization nis at 
T = 0. Based on unbiased QMC calculations for the 
three different dimerized models defined in Eq. ([T|) and 
Fig.lll we show that the curve TN{ms) exhibits a remark- 
able universality when properly normalized, not just close 
to the quantum-critical point but extending to strongly 
ordered systems. Our results give a parameter-free scal- 
ing function that can be compared with experiments. 

Quantum Monte Carlo calculations. — We have used 
the stochastic series expansioiifSSE) QMC method with 
very efficient loop updates [5|-|7| to calculate the squares 
(ml) and {m^^) of the z-components of the uniform and 
staggered magnetizations, 
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FIG. 2: (Color online) Procedures used to extract the critical 
temperature Tn- (a) and (b) show p^L and i?2, respectively, 
for the columnar dimer model at coupling ratio g = 3.444. 
The error bars are smaller than the symbols. Using polyno- 
mial fits to data for two lattice sizes, L and L -I- 2, crossing 
points between the curves are extracted. Results are shown 
in (c), along with fits of the form Tjv(L) = Tjv(oo) +a/ L™ (to 
the large-L data for which this form obtains). Extrapolations 
of the three quantities give Tn = 0.7996(3), 0.7996(6), and 
0.7999(5) for L — >■ oo, all consistent with each other within 
errors bars. 



where the phases 0; = ±1 correspond to the sublattices 
of the bipartite systems in Fig. [TJ The uniform suscepti- 
bility is X = {f^l)/{TN). We also study the Binder ratio, 
i?2 = ("^sz)/('^s2)^' 3,nd the spin stiffness constants p" 
in all lattice directions {a = x,y,z), p" = (PE{6a)/d6^, 
where E is the internal energy per spin and 9a a uniform 
twist angle imposed between spins in planes perpendicu- 
lar to the a axis. The stiffness constants can be related 
to winding number fluctuations in the simulations Q- 

We use standard finite-size scaling 0] to extract Tn- 
At Tn, the stiffness constants scale with the system 
length as p" oc L^"**, where the dimensionality d = 3. 
Thus, p"L should be size-independent at Tjv, while this 
quantity vanishes (diverges) for T > Tn [T < Tn)- In 
practice, this means that curves versus T (at fixed g) for 
two different system sizes L cross each other at a point 
which drifts (due to scaling corrections) toward Tn with 
increasing L. The dimensionless Binder ratio also has 



this kind of behavior and provides us with a different Tn 
estimate to check for consistency. Figs. [l^a,b) show ex- 
amples of these crossing behaviors for p^L and R2- The 
crossing points drift in different directions and bracket 
Tn- Fig.[2I^c) shows the L dependence of crossing points 
extracted from data for (L, L-f-2) system pairs, for R2 and 
two different stiffness constants. Power-law fits are used 
to extrapolate to infinite size. The mutual consistency of 
the Tn value so obtained using different quantities gives 
us confidence in the accuracy of this procedure. 

To extract the T = sublattice magnetization, we 
carry out simulations at temperature T = Ji/L. Note 
that, in a Neel phase with Tn > 0, any T{L) such that 
T{L 00) — )■ can be used for extrapolations to the 
thermodynamic limit and T = 0. Our choice is a natu- 
ral way to to scale the temperature since the lowest spin 
waves have energy oc 1/L. We also did some calculations 
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FIG. 3: (Color online) Extrapolation of the sublattice magne- 
tization obtained in simulations with T = Ji/L of the double- 
cube Heisenberg model at different coupling ratios g. The 
error bars are much smaller than the symbols. The fitting 
function used for L ^ oo extrapolations is a 
(where we exclude the linear term because it comes out very 
close to zero in fits including it). 

with T = 1/2L and obtained consistent extrapolated re- 
sults. Examples of the L dependence are shown in Fig. [3] 
for the double-cube model at several different coupling 
ratios. Taking into account rotational averaging in spin 
space, the final result for the sublattice magnetization is 
given by the L — >■ cxd extrapolated {ral^) (for which we 
use a polynomial fit, as shown in Fig.|3]); rris = y^3(mU). 

Universality of Tn versus mg. — Following the above 
procedures, we have calculated Tjv and rrig accurately 
for all three dimcr models at several coupling ratios g, 
from close to gc to deep inside the Neel phase. We graph 
T/v versus in Fig. HI is scaled by three differ- 
ent energy units; the intcr-dimcr coupling Ji in (a), the 
sum of couplings connected to each spin in (b), and 
the temperature T* at which the susceptibility exhibits 
a peak in (c). Before discussing these normalizations of 
T/v in detail, let us examine the reason for the linear be- 
havior, Tn oc TOs, seen in the QMC results for small [and 
in (b),(c) even quite large] m^. 

A semi-classical mean-field argument (inspired by the 
"renormalized classical" picture developed in two dimen- 
sions [H) leading to rUg is the following: To com- 
pute T/v in a classical system of spins of length S, one 
replaces the coupling of a spin So to the total spin of its 
neighbors i5, S^, by the thermal average JJ2s(^s)- 
In the presence of quantum fluctuations, this mean field 
seen by Sq is reduced, which is taken into account by a 
renormalization; (S^) — )■ {ms/S){Ss). The thermal fluc- 
tuations are, thus, added on top of the quantum fluctua- 
tions at r = 0, under the assumption that the quantum 
effects will not change appreciably for T > (i.e., the 
thermal fluctuations are assumed to be solely responsi- 
ble for further reducing the order). Note that Sq should 
not be renormalized here, but is computed as a thermal 
expectation value and should satisfy the self-consistency 
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FIG. 4: (Color online) The Neel temperature Tjv versus 
the sublattice magnetization for the three different dimerized 
models and with Tjv normalized in three different ways. Tjv 
is measured in units of (a) the inter-dimer coupling Ji, (b) 
the total coupling Js per spin, (c) the peak temperature T* 
of the susceptibility. A linear dependence obtains in all cases 
for small to moderate rUs, as indicate by fitted lines. Note 
that rris < 1/2 for S = 1/2. 

condition (S^) = (Sq)- The final magnetization curve is 
given by {ms/S){So). In this procedure of decoupling the 
classical and quantum fluctuations, one clearly effectively 
has J — > [rris/ S)J and, thus, T/v rrig. 

The assumption that the quantum renormalization fac- 
tor rris/ S is T-indcpendcnt up to Tjsi can be valid only 
if T/v is small. The energy scale in which to measure 
T/v when stating this condition should be dictated by the 
spin-wave velocity, which stays non-zero at the quantum- 
critical point [l7[ [i.e., not by the long-distance energy 
scale Ps{T = 0), which vanishes as g ^ gc and is un- 
related to the density of thermally excited spin waves]. 
A linear dependence is seen in Fig. 2] up to rather large 
values of rrig (where T/v ~ Ji)- A linear dependence was 
also recently found in the columnar dimcr model based 
on high-T expansions [l8| (with much larger error bars). 

Returning now to the issue of how to best normalize 
T/v, we note that in Fig. |4l[a), where the inter-dimer cou- 
pling Ji is used, the curve for the double-cube model is 
significantly above the other two. This is clearly because 
the constant Ji does not account for the different aver- 



4 




0.00 I 1 ' 1 ' 1 ' 1 ' 1 ' 1 ^ 

0.5 1.0 1.5 2.0 2.5 3.0 T 



3.0 


• Columnar Dimer 






Staggered Dimer 




2.5 


■ Double Cube 




2.0 






1.5 










(b) - 


1.0 







1 2 3 4 g 



FIG. 5: (Color online) (a) Susceptibility versus temperature 
of the staggered dimer model at different coupling ratios. The 
system size is L = 12, for which the peak height and location 
are already L — > oo converged, (b) The peak temperature 
versus the coupling ratio for the three different models. 

age couplings in the models. Using instead the sum of 
couplings connected to each spin, i.e., = 5 + g for the 
columnar and staggered dimers and 6 + 5 for the double 
cube (setting Ji = 1). the curves, shown in Fig. HJ^b), 
collapse almost on top of each other. Note that also the 
curves for the columnar and staggered dimers are closer 
to each other than in Fig. IH^a), although they have the 
same definition of J^. This can be the case because J, 
rescales the curves non- uniformly, since rus (g) and, there- 
fore, Js(ms), is different for the two models. The linearity 
of Tn/Js versus nis is also much clearer than before and 
extends all the way up to rris ~ 0.3. 

Although the data collapse is already quite good in 
Tn I Js , we can do even better when normalizing with a 
physical quantity that measures the effective lattice-scale 
energy. One such energy scale in antiferromagnets is the 
temperature at which the uniform magnetic susceptibil- 
ity X exhibits a peak. This peak is due to the cross-over 
from the high-T Curie form to the low-T weakly temper- 
ature dependent form typical of antiferromagnets. The 
peak temperature T*, thus, reflects the short-distance 
energy scale at which antiferromagnetic correlations be- 
come significant. T* is often used experimentally to ex- 
tract the value of the exchange constant, u sing , e.g., the 
"Bonner- Fisher" curve in one dimension [l9|. In spa- 
tially anisotropic systems such as the dimerized models 
we consider here, a natural assumption is that T* reflects 
an effective average coupling. In Fig. we show ex- 
amples of the susceptibility close to its peak, and in (b) 
we show the dependence of T* on g for all three mod- 
els. Normalizing Tfq with T* leads to remarkably good 
data collapse, as shown in Fig. SJc). Deviations from a 
common curve are barely detectable. Although we can- 



not prove that this function is really universal for all 3D 
networks of dimers, the results are very suggestive of this. 

Discussion. — The universal behavior implies that the 
T > disordering mechanism in the 3D Neel state is com- 
pletely governed by a single lattice-scale energy (which, 
as we have shown here, can be taken as the peak temper- 
ature T* of the susceptibility) and the T = sublattice 
magnetization TOs . The extended linear behavior seen in 
Figs. HlJbjC) shows that the quantum and classical fluc- 
tuations at T < Tn are completely decoupled all the 
way from g = gc (excluding gc itself, where Tn = 0) 
to quite far away from the quantum-critical point. De- 
pending on a lattice-scale energy instead of the quantum- 
critical spin stiffness, the linear behavior is not funda- 
mentally a quantum-critical effect. We have discussed 
the linearity and decoupling of the fluctuations in terms 
of a semi-classical mean- field theory, the validity of which 
implies that the quantum-critical regime 0] commences 
only above T/v. Deviations from linearity at larger TOs 
show that the quantum fluctuations are affected (become 
T-dependent) here, due to the high density of excited spin 
waves as r ^ T/v because is high. It is remarkable 
that this coupling of quantum and classical fluctuations 
also takes place in an, apparently, universal fashion for 
different systems. It would be interesting to explain this 
more quantitatively, by deriving the full function Tjq ver- 
sus nia analytically. Progress in the linear reg ime has 
been made recently in work parallel to ours [20[ . 

From a practical point of view, the data collapse of 
T/v /T* versus is very useful, because all the quantities 
involved can be measured experimentally and do not rely 
on microscopic details. The universal curve can be used 
to test the 3D Heisenberg scenario without adjustable 
parameters. The universality likely applies not only to 
dimer netwtorks, but also to systems where the quantum 
fluctuations are regulated in other ways. 
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